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We report an observation of the process J/ijj 7X(1835) —>■ 'yK^Kgr; at low KgKg mass with a 
statistical significance larger than 12.9o' using a data sample of 1.31 x 10® J/^|) events collected with 
the BESIII detector. In this region of phase space the K%K% system is dominantly produced through 
the/o(980). By performing a partial wave analysis, we determine the spin parity of the X(1835) to be 
J^^ = 0 The mass and width of the observed X(1835) are 1844 ± 9(stat)l25(syst) MeV/c® and 
192ti?(stat)tg (syst) MeV, respectively, which are consistent with the results obtained by BESIII 
in the channel J/?/) 'y'K'^'K~ri'. 


PACS numbers: 13.20.Gd, 13.66.Be, 14.40.Be 
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The non-Abelian property of quantum chromodynam¬ 
ics (QCD) permits the existence of bound states beyond 
conventional mesons and baryons, such as glueballs, hy¬ 
brid states, and multiquark states. The search for these 
unconventional states is one of the main interests in ex¬ 
perimental particle physics. One of the most promising 
candidates, the A(1835) resonance, was first observed in 
its decay to Tr~^TT~r]' in the process J/'0 —>■ by 

BESII [1]; this observation was subsequently confirmed 
by BESIII [2]. The discovery of the A(1835) has stimu¬ 
lated theoretical speculation concerning its nature. Pos¬ 
sible interpretations include a pp bound state [3] , a sec¬ 
ond radial excitation of the r]' [4], and a pseudoscalar 
glueball [5] . In addition, an enhancement in the invariant 
pp mass at threshold, X{pp), was first observed by BESII 
in the decay J/tp ^ jpp [6], and was later also seen by 
BESIII [7] and CLEO [8]. In a partial-wave analysis of 
J/'0 —>■ jpp, BESIII determined the of the X{pp) to 
be 0 ^ [9]. The mass of the X{pp) is consistent with the 
A(1835) mass measured in J/ip —>■ 77r+7r“?7' [2], but the 
width of the X{pp) is significantly narrower. 

To understand the nature of the A(1835), it is crucial 
to measure its and to search for new decay modes. 
Because of its similarity to J/ip —>■ yTr+Tr'r;', J/ip 
■jKKr] is a favorable channel to search for 
KKrj. In contrast to J/ip —>■ 'yK^K~r], there is no 
background contamination for J/ip —>■ from 

J/ip —>■ TT^KgKgT] and J/ip —>■ KgKgi/, which are forbid¬ 
den by exchange symmetry and CP conservation. There¬ 
fore, the channel J/ip —>■ jKgKgri provides a clean envi¬ 
ronment with minimal uncertainties due to background 
modeling. In this Letter, we report the first observa¬ 
tion and spin-parity determination of the X(I835) in 
J/ip ^ ’-fKgKg'r], where the Kg and r] are reconstructed 
from their decays to tt+tt” and 77, respectively. The 
analysis is based on a sample of (1310.6± 10.5) x 10® J/ij} 
events [10, 11] collected with the BESIII detector [12]. 

The BESIII detector is a magnetic spectrometer oper¬ 
ating at BEPCII, a double-ring e+e" collider with center 
of mass energies between 2.0 and 4.6 GeV. The cylin¬ 
drical core of the BESIII detector consists of a helium- 
based main drift chamber (MDC), a plastic scintilla¬ 
tor time-of-flight system, and a CsI(Tl) electromagnetic 
calorimeter (EMC) that are all enclosed in a supercon¬ 
ducting solenoidal magnet providing a 1.0 T (0.9 T in 
2012, for about 1087 x 10® collected J/ip) magnetic field. 
The solenoid is supported by an octagonal flux-return 
yoke with resistive plate counter muon identifier mod¬ 
ules interleaved with steel. The acceptance of charged 
particles and photons is 93% of the 47r solid angle, and 
the charged-particle momentum resolution at 1 GeV/c is 
0.5%. The EMC measures photon energies with a reso¬ 
lution of 2.5% (5%) at 1 GeV in the barrel (end caps). 
A GEANT4-based [13] Monte Carlo (MG) simulation soft¬ 
ware package is used to optimize the event selection cri¬ 
teria, estimate backgrounds, and determine the detection 
efficiency. 

Charged tracks are reconstructed using hits in the 


MDC. Because there are two with displaced vertices, 
the point of closest approach of each charged track to 
the e“''e“ interaction point is required to be within ±30 
cm in the beam direction and within 40 cm in the plane 
perpendicular to the beam direction. The polar angle 
between the direction of a charged track and the beam 
direction must satisfy [cos^J < 0.93. Photon candidates 
are selected from showers in the EMC with the energy 
deposited in the EMC barrel region (jcos^J < 0.8) and 
the EMC end caps region (0.86 < ] cost*] < 0.92) greater 
than 25 MeV and 50 MeV, respectively. EMC cluster 
timing requirements are used to suppress electronic noise 
and energy deposits unrelated to the event. 

Candidate J/ip —>■ ^KgKgr] events are required to 
have four charged tracks with zero net charge and at least 
three photon candidates. All charged tracks are recon¬ 
structed under the pion hypothesis. To reconstruct ATg 
candidates, the tracks of each 7r“®7r“ pair are fitted to a 
common vertex. candidates are required to satisfy 
I < 0.009 GeV/c^ and L/gl > 2, where 
L and ul are the distance between the common vertex 
of the 7r“'"7r“ pair and the primary vertex, and its error, 
respectively. The ■jjjKgKg candidates are subject to a 
kinematic fit with four constraints (4(7), ensuring energy 
and momentum conservation. Only candidates where the 
fit yields a X 4 c less than 40 are retained for further 
analysis. For events with more than three photon candi¬ 
dates, multiple J/ip ^ xKgKgT] candidates are possible. 
Only the combination yielding the smallest is re¬ 
tained for further analysis. Candidate J/tp ^ 'jKgKg-q 
events are required to have exactly one pair of photons 
within the ry mass window (0.51 < < 0.57 GeV/c^). 

Simulation studies show this criterion significantly re¬ 
duces the miscombination of photons from 3.20% to 
0.16%. The miscombination of pions is also studied and 
found to be negligible. To further suppress background 
events containing a tt®, events with any photon pair 
within a tt® mass window (0.10 < < 0.16 GeV/c^) 

are rejected. The decay J/ip (pKgKg with ())—>■ 777 
leads to the same final state as the investigated reaction 
J/tp ^ 'yKgKgrj. Therefore, events in the mass region 
\Mryn — m^\ <0.04 GeV/c^ are rejected. 

After applying the selection criteria discussed above, 
the invariant mass spectrum of KgKgr] shown in Fig. 1 
(a) is obtained. Besides a distinct rjc signal, a clear 
structure around 1.85 GeV/c^ is observed. The KgKg 
mass spectrum, shown in Fig. 1 (b), reveals a strong 
enhancement near the KgKg mass threshold, which is 
interpreted as the /o(980) by considering spin-parity 
and isospin conservation. The scatter plot of the in¬ 
variant mass of KgKg versus that of KgKg-q is shown 
in Fig. 1 (c). A clear accumulation of events is seen 
around the intersection of the /o(980) and the struc¬ 
ture around 1.85 GeV/c^. This indicates that the struc¬ 
ture around 1.85 GeV/c^ is strongly correlated with 
/o(980). By requiring <1.1 GeV/c^, the struc¬ 

ture around 1.85 GeV/c^ becomes much more prominent 
in the KgKgr] mass spectrum [Fig. 1 (d)]. In addition, 
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there is an excess of events around 1.6 GeV/c^. 



FIG. 1. Invariant mass distributions for selected events: In¬ 
variant mass spectra of (a) KgKgrj and (b) KgKg-, (c) scat¬ 
ter plot of versus (d) K%K%r] invariant 

mass spectrum for events with the requirement < 

1.1 GeV/c^. Dots with error bars are data; the shaded his¬ 
tograms are the non -77 backgrounds estimated by the rj side¬ 
band; the solid histograms are phase space MC events of 
J/'il) 'yKgKgrj with arbitrary normalization. 

Potential background processes are studied using a 
simulated sample of 1.2 x 10® J/ip decays, in which 
the decays with measured branching fractions are gen¬ 
erated by EvtGen [14] and the remaining J/il; decays 
are generated according to the lundcharm [15] model. 
Simulated events are subject to the same selection pro¬ 
cedure applied to data. No significant peaking back¬ 
ground sources have been identified in the invariant mass 
spectrum of KgKgrj. Dominant backgrounds stem from 
J/ij) iKgKgTT^ and J/ip —?> . These non- 

77 backgrounds are considered in the partial wave anal¬ 
ysis (PWA) by selecting events from data in the 77 side¬ 
band regions defined as 0.45 < < 0.48 GeV/c^ and 

0.60 < M-y-y < 0.63 GeV/c®, and they account for about 
2.5% of the total number of events in the 77 signal region. 

A PWA of events satisfying < 

2.8 GeV/c® and M^o^o < 1.1 GeV/c® is performed to 
determine the parameters of the structure around 1.85 
GeV/c^. These restrictions reduce complexities due to 
additional intermediate processes. The signal ampli¬ 
tudes are parameterized as sequential two-body decays, 
according to the isobar model: J/t/: —>• 7 A, X ^ Yr] 
or ZKg, where Y and Z represent the KgKg and 
Kgrj isobars, respectively. Parity conservation in the 
J/ilt ^ 'yKgKgri decay restricts the possible of the 
KgKgT] {X) system to be 0 1++, 2++,2 3++, etc. 


In this Letter, only spins J < 3 and possible 5'-wave or 
P-wave decays of the X are considered. The amplitudes 
are constructed using the covariant tensor formalism 
described in Ref. [16]. The relative magnitudes and 
phases of the partial wave amplitudes are determined 
by an unbinned maximum likelihood fit to data. The 
contribution of non-r; background events is accounted 
for in the fit by subtracting the negative log-likelihood 
(NLL) value obtained for events in the 77 sideband region 
from the NLL value obtained for events in the 77 signal 
region. The statistical signihcance of a contribution is 
estimated by the difference in NLL with and without 
the particular contribution, taking the change in degrees 
of freedom into account. 

Our initial PWA fits include an A'(1835) resonance in 
the /o ( 980)77 channel and a nonresonant component in 
one of the possible decay channels /o(980)77, /o(1500)77 or 
/2 ( 1525 ) 77 . All possible combinations of the A(1835) 

and the nonresonant component are tried. We then ex¬ 
tend the fits by including an additional resonance at lower 
KgKgr] mass. This additional component, denoted here 
as the A (1560), improves the fit quality when it is al¬ 
lowed to interfere with the A(1835). Our final fits show 
that the data can be best described with three compo¬ 
nents: A(1835) ^ /o(980)77, A(1560) /o( 980)77, and 

a nonresonant /o(1500)77 component. The of the 
A(1835), the A(1560), and the nonresonant component 
are all found to be O The A(1835), A(1560), and 
/o(1500) are described by nonrelativistic Breit-Wigner 
functions, where the intrinsic widths are not energy de¬ 
pendent. The masses and widths of the A(1835) and 
A(1560) are derived by scanning each over a certain 
range. The /o(1500) mass and width are fixed to the 
values reported in Ref. [17]. The /q( 980) is parame¬ 
terized by the Platte formula [18], with the parameters 
fixed to the values reported by BESII [19] in the channels 
J/ 7/7 —>■ ^7r+7r“ and J/ip ^ cpK'^K~ . The scan returns a 
mass and width of the A(1835) of 1844 ± 9MeV/c^ and 
MeV, respectively. The mass and width of the 
A(1560) are determined to be 1565 ± 8MeV/c® and 
4511^ MeV, respectively. Using a detection efficiency 
of 5.5%, obtained by a MG sample weighted by par¬ 
tial wave amplitudes, the product branching fraction of 
J/iP ^ 7A(1835) and A(1835) ^ ^ 0^077 (S^dsas)) 
is calculated to be (3.311o 3 q) x 10“®, where the decay 
A(1835) KgKgr] is dominated by /o(980) production. 
The KgKgT], KgKg, Kg-q mass spectra and the distri¬ 
butions of the J/ip, KgKgq and KgKg decay angles are 
shown in Fig. 2. Overlaid on the data are the PWA 
fit projections, as well as the individual contributions 
from each component. The x^/’^bin value is displayed 
on each figure to demonstrate the goodness of fit. We 
evaluate the significance by applying the likelihood ratio 
test, performing a separate fit for every systematic vari¬ 
ation detailed below. The most conservative statistical 
significances of the A(1835) and A(1560) are 12.9(7 and 
8.9(7, respectively. 

Various fits are performed by changing the and 
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FIG. 2. Comparisons between data and PWA fit projections, (a), (b), and (c) are the invariant mass distributions of K%K%rj, 
KgKg, and KgT] (two entries/event), respectively, (d)-(f) are the angular distributions of cos6, where 6 is the polar angle of (d) 
7 in the J/ijj rest system; (e) 77 in the K%K%ri rest system; and (f) K% in the K%K% rest system (two entries/event). The dots 
with error bars are data, the solid histograms are the PWA total projections, the shaded histograms are the non -77 backgrounds 
estimated by the rj sideband, and the short-dashed, dash-dotted, and long-dashed histograms show the contributions of A(1835), 
X(1560), and the nonresonant component, respectively. 


decay mode of the nonresonant component compared to 
the nominal solution described above. The NLL value 
of a fit with a I’*"’'' nonresonant /o(1500)r7 componeirt 
is only worse by 0.8 compared to the nominal solution, 
which indicates that we canirot distiirguish between the 
two spin assignirrents of the nonresoirairt component with 
our present statistics. This ambiguity introduces large 
systematic uncertainties in the Sx(r 835)5 since the inter¬ 
ference between the Al(1835), Ar(1560), and the nonres¬ 
onant component depends on the spin assignment of the 
latter. To establish the of the X(1835), we per- 
form a series of PWA fits assuming alternative hy¬ 
potheses for both the A(1835) and the nonresoirant con¬ 
tribution. For the iroirresonant contribution, we also test 
several possible decay channels [/o(980)77, /o(1500)77, and 
/2 ( 1525 ) 77 ] in turn. For each nonresonant component as¬ 
sumption, the A(1835) 0 ^ hypothesis is significarrtly 
better than the l+“'' or 2 ^ hypotheses, with the NLL 
value improving by at least 41.6 units. Analogously, we 
perform the same series of PWA fits for the A(1560). 
Again the 0 ^ hypothesis for the A(1560) always yields 


a significantly better fit result than other assign¬ 

ments, with the NLL value improving by at least 12.8 
units. 

We evaluate the contributions from additional well- 
known resonances by adding them individually to the 
fit. We consider all possible combinations for X 
and its subsequent decay products Y and Z as given 
in Ref. [17]: for X, this includes ?7(1760), 77(2225), 
/i(1510), 772(1870), /2(1810), /2(1910), /2(1950), 

/2(2010), /2(2150), /2(2300), /2(2340), /j(2220); for T, 
/o(980), /o(1500), /o(1710), /2(1270) and /2(1525); for 
z, a:*(14io), a:*(1680), a:^(1430), a:^(i950), a: 2*(1430), 
and Arj(1980). Additional nonresonant contributions 
with various and decay modes are studied as well. 
The statistical significances of the additional contribu¬ 
tions are smaller than 5 (t. In order to check the pos¬ 
sible contribution from a nonresonant KgKg process, 
we add a A(1835) {K^K^)s'rj process into the nom¬ 
inal solution, where {KgKg)s refers to a nonresonant 
KgKg contribution in a relative S wave. We find the 
resulting significance of the A(1835) — fo{980)ri pro- 
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cess and the X(1835) —)• {KgKg)sr] process to be 6.8tT 
and 1.6 ct, respectively, so we do not include the latter 
process in the nominal solution. We also test a fit by 
changing the decay mode of the X(1560) in the nom¬ 
inal solution from fo{980)r] to {KgKg)s'r]] the fit with 
X(1560) — {KgKg)sr] has almost the same fit quality as 
the nominal solution. Therefore, with the present statis¬ 
tics, we cannot draw a conclusion about the X(1560) de¬ 
cay mode. The largest differences in masses and widths 
of the X(1835) and X(1560) and the product branching 
fraction i3x(i835) between all above alternative fits and 
the nominal solution are taken as systematic uncertain¬ 
ties from the components in the nominal solution. 

For the measurements of the masses and widths of the 
X(1835) and X(1560) and the product branching frac¬ 
tion Sjf(i 835 )i we include the following sources of sys¬ 
tematic uncertainties in addition to the sources discussed 
above: we change the KgKg mass range to < 

1.05, 1.15 and 1.20 GeV/c^; we change the /o(980) 
mass and coupling constants in the Flatte formula to 
other experimental measurements [ 20 - 22 ] ; we change the 
/o(1500) mass and width by one standard deviation [17]; 
we increase and decrease the non-r; background level by 
one standard deviation; we change the parameterization 
of the X (1835) and X(1560) line shape to a Breit-Wigner 
function whose intrinsic width is energy-dependent [23]; 
and we replace the X(1560) by ?7(1405) or ?7(1475). For 
the systematic errors of the product branching fraction 
6 x(i 835 )i we also consider the following additional un¬ 
certainties. The Kg reconstruction efficiency is stud¬ 
ied using two control samples of J/ip —>• K*^K^ and 
J/^p (j)K°gK^ TT^, while the photon detection efficiency 
is investigated based on a clean sample aiJ/ip^pir. The 
differences between data and MC simulation are 1.0% 
for each Kg and 1.0% for each photon [24]. A control 
sample of J/'ij} —>■ jKgKgTT^ is selected to estimate the 
uncertainty associated with the 4(7 kinematic fit. The ef¬ 
ficiency is the ratio of the signal yields with and without 
the kinematic fit requirement xlc ^ '^ 0 . The difference 
between data and MC simulation, 1.5%, is assigned as 
the systematic uncertainty. We also consider the uncer¬ 
tainties from the number of J/ip events [10, 11] and the 
branching fractions of K^ tt+tt” and rj —>■ 77 [17]. We 
change the mass and width of Ar(1835) or 7f(1560) by 
1 standard deviation of the statistical uncertainty. The 
individual uncertainties are assumed to be independent 
and are added in quadrature to obtain the total system¬ 
atic uncertainties as presented in the Supplemental Ma¬ 
terial [25]. 

In summary, a PWA of J/'0 —)■ jKgKgr] has been 
performed in the mass range < 2.8 GeV/c^ af¬ 
ter requiring < 1-1 GeV/c^. The PWA fit re¬ 

quires a contribution from Ar(1835) —>■ KgKgr] with 
a statistical significance greater than 12.9cr, where the 
7f(1835) — KgKgT] is dominated by /o(980) produc¬ 


tion. The spin parity of the X(1835) is determined 
to be 0 The mass and width of the Ai(1835) 
are measured to be 1844 ± 9(stat)l25(syst) MeV/c^ and 
192ti7(stat)l43(syst) MeV, respectively. The corre¬ 
sponding product branching fraction 6 x(i 835 ) is mea¬ 
sured to be ( 3 . 3 lto; 3 o(stat)l 4 ' 2 g(syst)) x 10“^. The mass 
and width of the Ai(1835) are consistent with the values 
obtained from the decay J/ip ^ yTr+Tr^ry' by BESIII [2]. 
These results are all first-time measurements and provide 
important information to further understand the nature 
of the X(1835). 

Another 0“+ state, the Ar(1560), also is observed 
in data with a statistical significance larger than 8.9tT 
and is seen to interfere with the 7f(1835). The mass 
and width of the Al(1560) are determined to be 1565 ± 
8 (stat)lg 3 (syst) MeV/c^ and 45ti3(stat)t28(syst) MeV, 
respectively. The mass and width of the X(1560) are 
consistent with those of the ?7(1405) and 7(1475) as given 
in Ref. [17] within 2.0a and 1.4 ct, respectively. Present 
statistics do not allow us to conclusively determine if the 
7f(1560) is the same state as the 7(1405)/7(1475) or a 
new meson. More statistics in this analysis and an am¬ 
plitude analysis of J/ip 777 r° 7 r° and J/ip —>■ yATgAT^Tr^ 
processes may help to understand the nature of the 
X(1560). 
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